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Abstract: We study the gauge/gravity duality for supersymmetric SU{N) Yang-Mills 
theory in l-|-0 dimension with sixteen supercharges using lattice simulations. The con¬ 
jectured duality states that the gravity side is described by N DO-branes in type IIA 
superstring at large N, and the thermal gauge theory reproduces the black hole thermo¬ 
dynamics at low temperature. In this paper, we explain the Sugino lattice action used 
in the simulations in detail, and examine the supersymmetric Ward-Takahashi identity to 
confirm the restoration of supersymmetry in the continuum limit. We also estimate the 
internal energy of the black hole from the lattice results for N = 14, 32, and find that it 
smoothly approaches the prediction of the gravity side as the temperature decreases. 
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1 Introduction 

The gauge/gravity duality is an equivalence between strongly coupled gauge theory and 
the theory of gravity on a curved space, which is regarded as a generalization of the original 
AdS/CFT correspondence [1]. At present, it is widely accepted, and many related studies 
have been carried out in various fields of physics [2]. However, it is not a theorem but 
a conjecture. In the verihcation of the duality conjecture, numerical simulations of the 
lattice gauge theories will play a key role. In this paper, we perform lattice simulations 
of one dimensional supersymmetric SU{N) Yang Mills theory with sixteen supercharges 
(’’the BFSS model” [3]) and examine the duality in this system. 

The gravity dual of this theory is described by N DO-branes in type IIA superstring 
in the large N limit [4], and the behavior of the thermal gauge theory at low temperature 
reproduces that of a black hole [5]. However, at low temperature, the gauge theory is a 
strongly coupled one in which known analytic techniques such as the perturbation theory 
become less effective. In addition, localization methods [6] and other analytic treatments 
[7] for BPS states are not applicable to the present case because supersymmetry is broken 
by temperature. Therefore, some numerical approaches are needed to examine the duality. 

Numerical simulations of this gauge theory have been carried out independently by two 
collaborations: M. Hanada, J. Nishimura and others [8], [9], [10], [11] and S. Catterall and 
T. Wiseman [12], [13], [14]. ^ Some regularization schemes are needed to put the theory 

^ Lattice simulations of two dimensional SYM with sixteen supercharges have been performed in [15], 
[16]. 
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on a computer. Hanada et al. used the momentum sharp cutoff, while Catterall-Wiseman 
employed the standard lattice theory. In both the theories, all of the supersymmetries 
are broken not only by temperature but also by the regularization schemes. However, as 
they expected, there would be no serious problems because the SUSY breaking effects of 
the regularization schemes disappear in the continuum limit thanks to the UV finiteness 
of one dimensional gauge theory. Some obtained results imply the validity of the duality. 
Nevertheless, further simulations are required to verify it with high accuracy. 

In this paper, we aim to verify the duality conjecture from lattice simulations of the 
gauge theory side. In particular, we employ the Sugino lattice action with keeping two 
supercharges on the lattice [17], [18], [19]. In section 2, we briefly review the continuum 
theory and its twisted action, after which the lattice action is defined from the twisted 
one. Section 3 describes some algorithms and techniques used in the simulations. We show 
numerical results of the supersymmetric Ward-Takahashi identity in section 4, then present 
the internal energy of the black hole in section 5. Finally, section 6 is devoted to a summary. 


2 One dimensional SYM with 16 supercharges 


In this section, we introduce one dimensional supersymmetric Yang-Mills theory with six¬ 
teen supercharges (the BFSS model [3]). Using two nilpotent supercharges Q±^ the contin¬ 
uum action is expressed as a Q±-exact form, and the Sugino lattice action is defined from 
the exact one [19]. ^ 


2.1 Continuum theory 

Supersymmetric Yang-Mills theory in l-|-0 dimension with sixteen supercharges contains 
a gauge field Aq, nine scalars Yj (i = 1, • • • ,9) and sixteen fermions '0^ (a = 1, • • • , 16). 
The continuum euclidean action is given by 


Srnn.t. — 


^ jdt tr|i(17oY,)2-i[Y„Y,f 




( 2 . 1 ) 


where A is the ’t Hooft coupling constant and 7 * are real symmetric matrices satisfying 
the nine dimensional Clifford algebra, { 7 i, 7 j} = 25ij. Also, the covariant derivative Dq is 
defined as Dq ■ = dt ■ +i[Ao, • ]. 

In one dimension, the gauge held has no kinetic terms and appears only through the 
covariant derivative. However, as well as four dimensional gauge theory, it leads to an 

^ In this paper, the gauge group is SU{N) and, unless otherwise noted, all of the fields are expressed 
as the matrix-valued functions. The gauge group generators r“(a = 1, • • • — 1) are traceless hermitian 

matrices which satisfy the Lie algebra, [r“,r^] = ifabcT"^, with the normalization tr(r“r*’) = Sab, The 
scalar and fermi fields, as well as the gauge field To, are expanded as (p(t) = Moreover, in 

this section, we impose the periodic boundary condition on all of the fields to study the SUSY invariance, 
while in the following sections the fermions satisfy anti-periodic boundary condition to study the theory at 
finite temperature. 
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invariance of the theory under the gauge transformations, 

SujAo = -Dqu, ( 2 . 2 ) 

for G and u is an infinitesimal gauge parameter. In this case, one can fix the 

gauge completely: Aq = 0 for no boundary conditions, or for periodic time where 

tc is a constant corresponding to the Polyakov line. Nevertheless, we define the theory 
without any gauge fixing to treat both SUSY and gauge invariance manifestly. 

The supersymmetry transformations are given by 

<^ 5^0 ~ ^a'4’at 

(^-3) 

1 % 

~ 2 (^b)a/3 ~ 

where are sixteen fermionic parameters. We can show that the action (2.1) is invariant 
under the transformations (2.3) by using the Fierz identity for the fermions and the Jacobi 
identity for the scalars. The well-known Noether’s theorem tells us that the corresponding 
supercharges Qa are the generators of the infinitesimal transformations, = ^aQa- 

This one dimensional theory can be obtained by a dimensional reduction of four di¬ 
mensional AA = 4 SYM, or its topologically twisted theory [20]. In the topological field 
theory, Lorentz group is twisted by the internal group of AA = 4 SYM. The transformation 
laws of the fields and the supercharges under the twisted Lorentz group are different from 
those under the original one. The two linear combinations of the sixteen supercharges Q± 
are scalars which are nilpotent up to the gauge transformations. Furthermore, one can 
express the action with a topological term as a (5±-exact form. The Q^-invariance of the 
action follows from the nilpotency of Q±. 

In order to define the Q^-exact action directly from (2.1), we only have to change the 
field variables to those of the twisted theory. We use Bi, C, 4>± and 'i/’i/ii X±i as the 
bosonic and fermionic twisted fields, respectively, for /i = 0,1, 2, 3 and i = 0,1, 2. Under the 
four dimensional twisted Lorentz transformations, A^,'ip±^ transform as vectors, C,(j)±,r]± 
remain unchanged, (namely are scalars), and the other bosons and fermions transform as 
self-dual tensors. Bi,x±i are three independent components of the bosonic and fermionic 
tensors, respectively. In one dimension, the indices /r, i should be regarded as the names 
of the helds. 
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The bosonic twisted fields A^,Bi, C, 4>± are given by 


X/i — [fj, — 1 , 2 , 3 ) 

^4 = -B2, 

^5 =^1, 

Xe = -Bo, (2.4) 

= Ic, 

^9=2 ('^+ + 'z^-) 5 

with Aq unchanged. A^,Bi,C are hermitian and (</>+)^ = —<p-- 

The fermionic twisted variables V’lt/i) X±i and the original variables ^jja are related 
to each other by 




^ 0-0 + f ??+ ^ 

'4’2 


0-1 - (X+2 

V’3 


0-2 + *X+1 

■04 


0-3 - «X+0 

05 


0+0 + f ??- 

06 


0+1 + *X-2 

07 


0+2 - *X-1 

08 

1 

0+3 + *X-0 

09 

~V2 

-i(0_o - §r/+) 

010 


-i(0_i + iX+2) 

011 


-i( 0_2 - fx+i) 

012 


-i(0_3 + iX+o) 

013 


-00+0 - §??-) 

014 


-00+1 - fx_2) 

015 


-00+2 + *X-l) 

V016/ 


\-00+3 - *X-0)/ 


with the appropriate gamma matrices. ^ 

® We take the following gamma matrices, 

71 = 0-2 (g) 1 ® (73 (g) (72, 74 = (73 ® (72 (g 1 ® (72, 77 = (73 (g (7l (g 1 (g 1, 

72 = (72 (g 1 (g (72 (g 1, 75 = (TS <g (^2 ® (72 (g (73, 78 = (^3 <g (^3 ® 1 <g 1, (2.6) 

73 = 0'2 g 1 0 C 7 l g (72 , 76 = ((3 g (^2 g (72 g (71, 79 = (7l g 1 g 1 g 1, 

where ai are the Pauli matrices. The gamma matrices act on the fermions as 

(71 )cKl CK2 ^((2 (((2)7*2 ® ^1112 g g (((2) 7 (2 "^(*^2 (^T) 

for a = 8(i — 1) + 4(j — 1) + 2{k — 1) + Z. 
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Let us define two supercharges Q± as follows: 


Q+ 

Q- 


“ 1 “ 


V2 


{Qi + iQg). 


From (2.3), we can read Q^-transformation laws, 



Q±'ilj±ii = -iDf,4>±, 

Q±Bi 

Q±x±i = [Bi,(t)±], 

Q±C = r]±, 

Q±v± = [C,4>±], 


Q±Vt ~ ['A=F> ‘Ai]) 

Q±4>± = 0, 

Q±X=f* = ± Hi 

± Hf^, 


( 2 . 8 ) 

(2.9) 


( 2 . 10 ) 


Q±Hi = ±{[x^i,(l)±] + ^[x±i,C] + ^[Bi,r]±]), 

Q±H^ = ± 4>±] + C] - ^D^r]±) , 

where Dq is the covariant derivative and D^{iJ, = 1, 2, 3) are the commutators •]. 

Note that are used for /r = 0,1, 2,3 because the twisted theory is originally defined in 
four dimensions. Moreover, seven auxiliary fields Hi and are introduced to make the 
transformations satisfy 


Q± = > {Q+’ Q-} = -i^c, 


( 2 . 11 ) 


where the right-hand sides are the gauge transformations (2.2) with the field-dependent 
gauge parameters. 

Performing the change of variables (2.4) and (2.5) with the gamma matrices (2.6), ® 
the action (2.1) can be expressed as the Q±-exact form, 


Scant. — dt tr ^ 2,iBi -I- 2 ^(-ijkBi\BjB}^ 


-V’+mV’-m - X+iX-i - jV+V 


(2.13) 


"^The fields are imaginary, in general. An imaginary function 0 can be expressed as cf) = (fjR + icj)i 
where (j)ji and (j)i are real functions. Similarly, the gauge transformation means = 5^^ + iS^j. 

® In addition, we added the following Gaussian integrals to the action, 


S'-S-S'-l-y J dttrj - i ^Fis + ^SijkFjk^ - ^eijk[Bj, Bk]J 

+ (^iii — i{tijkDjBk + DsBi)^ + -|- iDiBi^ . 


( 2 . 12 ) 


- 5 - 



where 


Foi = DqAi, Fij = i[Ai,Aj], {i,j = 1,2,3). (2.14) 

These F^i, are originally the field tensors in four dimensions and become the covariant 
derivatives and the commutators after the dimensional reduction. 

From the nilpotency (2.11), the action (2.13) is trivially invariant under the Q±- 
transformations, without a usage of the Leibniz rule. The other fourteen supersymme¬ 
tries and the global SO{9) symmetry are also exact symmetries of (2.13) because we only 
performed the change of variables from the original variables to the twisted ones. 

2.2 Lattice theory 

We construct the theory on a finite lattice of size L with periodic boundary condition. The 
sites of the lattice are labeled by the integers t = 1, • • • ,L and the lattice spacing o is set 
to unity without loss of generality. The adjoint scalars and fermions live on the sites, while 
a lattice gauge field U{t) G SU{N) lives on the links. 

The lattice gauge transformations are given by 

d^U{t) =-iVouj{t)U{t), =-i[(p{t),uj{t)], (2.15) 

where ip represents all of the scalars and fermions, and uj is an infinitesimal gauge parameter 
defined on the sites. Under the gauge transformations, the covariant forward difference 
operator 


VQipit) = U+ 1)U ^{t)-(p{t) 


(2.16) 


transforms as (p{t) itself. 

By replacing the integral and the covariant derivative in (2.13) and (2.14) with the 
summation over the sites and the covariant forward difference operator (2.16), respectively, 
we can introduce the following lattice action, ® 


( 2 . 

where 


N V—\ f / 1 \ 1 

S = Q+Q_— ^ tr |-2fBi ^Fis + ij^f^ijkFjkj - -eijkBi[Bj,Bk] 


-'i/’+M'i/’-M - X+iX-i - tWW 


Foi = Vo^i, Fij = i[Ai,Aj\, = 1,2,3). (2.18) 

The ’t Hooft coupling A has mass dimension three, and therefore Aq = Ao^ is the dimen¬ 
sionless ’t Hooft coupling constant. The continuum limit of this theory is obtained by 
letting Ao to zero with hxed A. 

® This is a one dimensional version of the action given in [19], which is therefore referred to as the Sugino 
lattice action in this paper. 
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The lattice Q±-transformations 


Q±U = i'ilj±oU, 

Q±Af, = Q±^P±im = -iV/,(/>± + i5;,oV'±oV'±o, 

Xiti) QitXiti 

Q±c = r]±, Q±^l± = [C,4>±], 

Q±^t ~ Q±^=f ~ [4'^t4‘±]^ 

Q±(i>± = 0, 


Q±XTi = iP, Bi] ± Hi, 

Q±4>Tf^ = + i5^o{V'+o, V^-o}, 



(2.19) 


are the same with the continuum ones with the exception of the fields with /i = 0. The 
terms proportional to 6^q are higher order corrections. These Q± satisfy 


Q% = i5^^, {Q+,Q-} = -i6c 


( 2 . 20 ) 


even on the lattice [19]. There is, however, nothing surprising about this if we recall the 
BRST transformation. In the BRST transformation, once the transformation law of a low 
dimensional field is given, then that of the higher dimensional field is uniquely determined 
from the nilpotency. Similarly, even on the lattice, once the Q±-transformation laws of the 
low dimensional helds {Q±U, Q±Bi, etc.) are given, then the others {Q±'4’, Q±X±j etc.) 
are uniquely determined from (2.20). 

The lattice action (2.17) possesses the (5j--invariance and the gauge invariance. On 
the other hand, the other fourteen supersymmetries and the S'©(9) global symmetry are 
broken at a finite lattice spacing. Instead of the SO{9) symmetry, the lattice theory has, 
at least, a SU{2) global symmetry corresponding to the interchange of Q+ and Q-. In 
the naive continuum limit, the lattice (5±-transformations (2.19) reproduce the continuum 
ones (2.10) and the lattice action reproduces the correct continuum action (2.1) via (2.13). 

After performing Q± in (2.17), we find that the lattice action has a four-fermi inter¬ 
action. 



( 2 . 21 ) 


which is of the order of the cut-off. The four-fermi term is not suited for the numerical 


simulations. In order to express the four fermion interaction as a fermion bilinear form, as 
in the NJL model, we introduce an auxiliary field a{t) and replace (2.21) with 



( 2 . 22 ) 


In fact, this reproduces (2.21) by integrating a. 
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3 Simulation details 


We used the standard rational Hybrid Monte Carlo method for the numerical simulations 
[21], As explained in section 2.2, the auxiliary held a was introduced to express the fermi 
action as a fermion bilinear from. Therefore, the HMC generated the conhgurations of the 
eleven bosonic variables, C/, Ai, • • • , Ag and o. Moreover, as explained below, we prepared 
the fermion pfaffian by using the pseudo fermion method with a rational approximation. 
For the phase of the pfaffian, we used the phase reweighting method for the SUSY WTI in 
section 4 and quenched it for the internal energy of the black hole in section 5. ^ 

The dynamical effects of the fermions can be included in the simulations through the 
pseudo fermion method. In the present model, the integration of the fermions yields the 
pfaffian pf(D) which is complex values in general. We treated the absolute value and the 
complex phase of the pfaffian, individually. Since |pf(i4)| = ^ the absolute 

value of the pfaffian can be given by 


|pf(F))| = y T)</)tT)^exp|-(/.^(T)tT»)-3^} , (3.1) 

where </> is a complex pseudo fermion, and the 4th-root of can be approximately 

prepared by a rational expansion. 




M 

ao + Y^ 

i=l 


CXi 

DWTWi' 


(3.2) 


where the degree M and the coefficients Oj, /3j depend on the range of the eigenvalues of 
and the accuracy of the approximation. We examined the range of the eigenvalues at the 
thermalization step in each simulation parameter and chose M, Oi , [3i of the approximation 
within the relative errors of 0(10“^^) [22]. Then, we fixed them during the configuration 
generation. We justified our initial guess for M,ai,l3i by computing the maximum and 
minimum eigenvalues in each trajectory. In addition, we computed the inversions of D 
with the shifts jSi in (3.2) by using the multi-mass CG solver [23]. The effect of the phase 
of the pfaffian was included in the results of the SUSY WTI using the phase reweighting 
method. Meanwhile, we merely ignored it when we measured the internal energy of the 
black hole. 

As already reported in [8], [13], [24], the HMC achieved the thermalization at high 
temperature, while at low temperature it was unstable and we observed that the magnitudes 
of the scalar fields have been monotonically increasing against the trajectories. Once such 

^ This is purely due to some technical reasons. The numerical calculation of the pfafRan phase is 
demanding and costs 0{K^) for K x K matrices. We have to use sufficiently large N for testing the gauge 
gravity duality in the large N limit and avoiding the instability related to the flat directions. If we use a 
lattice of size L = 16 and N = 16, then the size of the Dirac operator becomes K = {N'^ — 1) x 16 x L ~ 
65, 000, which is too large for the calculation of the phase. Thus, we did not compute the pfaffian phase 
when we measured the internal energy of the black hole. On the other hand, = 2, 3,4 are sufficient to 
confirm the restoration of SUSY from the SUSY WTI. For such small N, we computed the phase every 
configuration. 
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phenomena occurred, the thermalization was not achieved. ® Therefore, in the parameter 
regions with the instability, we added a mass term for the scalar fields. 


'mass 



(3.3) 


to the action, where /io(= is the dimensionless mass. This term was used to examine the 
SUSY restoration as presented in the next section. Meanwhile, in section 5, we measured 
the internal energy of the black hole without the mass term by taking sufficiently large N 
to avoid the instability. 

The trajectory length was r = 0.5, and the number of step of the molecular dynamics 
was chosen such that the acceptance rate was between 80% and 90%. After discarding the 
first O(IO^) trajectories for the thermalization, we stored 1 configuration every 10 trajec¬ 
tories. 

4 Restoration of supersymmetry 

In this section, we see that the effect of the lattice spacing vanishes and supersymmetry 
is restored in the continuum limit by numerically examining the supersymmetric Ward- 
Takahashi identity with the SUSY breaking term (partially conserved super current). 

The present lattice action has two exact supercharges, but the remaining fourteen 
supercharges are broken by the lattice cut-off. In the classical continuum limit, the broken 
supersymmetries are trivially restored, while in the quantum theory the restoration does 
not occur in general because some SUSY breaking relevant operators can be generated 
by the UV-divergences. Fortunately, in one dimensional gauge theories, the perturbative 
power-counting tells us that no SUSY breaking relevant operators exist, in other words, the 
theories are UV-finite. Hence, the supersymmetries are restored in the quantum continuum 
limit. However, it is not clear whether this argument surely holds beyond the perturbation 
theory and at finite temperature. Therefore, we need to show that the restoration does 
occur. 

4.1 Supersymmetric Ward-Takahashi identity 

In the following, we derive the supersymmetric Ward-Takahashi identity in the continuum 
theory at finite temperature using the path integral formulation. 

We define the partition function and the expectation value of an operator O as follows: 



(4.1) 


(4.2) 


“ The threshold temperature, where we could observe such phenomena, depended on N, and the larger 
N were more stable. This instability comes from the flat directions of the boson action. 
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where (p represents the whole field, and the continuum action is given by 

S'cont. = Scant. + ^ [ dt Y^tr{Xf{t)} . (4.3) 

d i=l 

We added the mass term to the action. Hereafter, we take a fermionic one point operator 
Op{s) as O. 

For the integration variables, we make a change of variables, which is given by the 
localized supersymmetry transformations, (2.3) with ^ —>• i{t). Note that the local pa¬ 
rameter ^(t) satisfies the anti-periodic boundary condition. Under the change of variables, 
Scant, changes to the total derivative of the supercurrent, 

N 

da = ^ [(7i)a/3 tr (ip^DoXi) - i(Sjj)«/3 tr{'ilJi3[Xi,Xj])] , (4.4) 

and the mass term yields a supersymmetry breaking term, 

N 

Xa = y(7i)a/5 trjWV’/?} ■ (4.5) 

The integral measure is invariant under the change. Thus, we Hnd 

/ ida{t)Op{s)) = J dt ^a{t) {Ya{t)Op{s)) 

- J dt6{t - s)Ca{t) {QaOpis)), (4.6) 

where the supercharges Qa in the last term give the infinitesimal supersymmetry transfor¬ 
mations. Since ^ in (4.6) is arbitrary, one can obtain the SUSY Ward-Takahashi identity, 

^ {Ja{t)Op{s)) = pt^ {Ya{t)Op{s)) - 6{t - s) {QaOpis)) (4.7) 

This derivation of the SUSY WTI is independent to the boundary condition, hence we 
find that (4.7) holds for both periodic and anti-periodic boundary conditions. From a 
lattice version of the supersymmetric Ward-Takahashi identity, one can investigate the 
supersymmetry breaking effect of the lattice spacing, as distinguished from those of the 
temperature and the mass [24]. 

4.2 Numerical results of SUSY WTI 

In this section, we present the numerical results of the SUSY WTI, including the effect 
of the complex phase of the pfaffian, for SU{2) and at the temperature = 1. The 

parameters set used in the simulations is presented in table 1. At this temperature, we 
observed the instability related to the flat directions. Therefore, as explained in section 
3, we added the mass term to the lattice action and measured the SUSY WTI with the 
breaking term Y. For the simplicity of the explanation, we give the results in the physical 
unit (A = 1). 
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SU(2) 

T/X^/^ = 

1 

L 

a 


traj. 

8 

0.125 

0.01 

15,000 

8 

0.125 

0.02 

15,000 

8 

0.125 

0.05 

15,000 

12 

0.0833 

0.01 

15,000 

12 

0.0833 

0.02 

15,000 

12 

0.0833 

0.05 

15,000 

16 

0.0625 

0.01 

15,000 

16 

0.0625 

0.02 

15,000 

16 

0.0625 

0.05 

15,000 


Table 1. Simulation parameters and the HMC trajectories. The lattice spacing a = (TL) ^ and 
the mass are presented in the physical unit (A = 1). 


As a lattice counterpart of the SUSY WTI, we chose Y (4.5) as the source operator O 
and calculated the following sixteen ratios, ® 


{dJa{t)Y^{0)) 
(Y„(t)y„(0)) ’ 


for a = 1, • • • , 16, 


(4.9) 


where d is the symmetric difference operator, dJa{t) = [Jait + 1) — Jo(t — l)]/2, and the 
lattice supercurrents Ja are prepared by replacing the covariant derivative in (4.4) with 
the covariant forward difference operator. These sixteen ratios correspond to the sixteen 
supercharges. Note that there are two exact supercurrents corresponding to Q±. Since their 
concrete forms are complicated, we used this simple definition for all of the supercharges. 
The calculation of the inversion of the Dirac operator is needed to estimate the ratios. We 
used LAPACK [25] to compute the inversion. Averaging the ratios over the lattice sites, 
we reduced the statistical errors. 

Figure 1 shows the ratios against time on which the supercurrents are defined for 
L = 16 and /i^ = 0.01. The corresponding lattice spacing is a = 0.0625. There are sixteen 
points corresponding to the sixteen supercharges on the same time-slice. Clear plateaus 
are seen for 3 < t < 13. All of the sixteen plateau values are = 0.01 within the errors. 
In other words, the SUSY WTI holds near the continuum limit. 

On the other hand, near t = 0 (and L), the ratios deviate from the plateaus because 
of the contact term in (4.7). The contact terms are smeared on the lattice, therefore the 
ratios near t = 0 (even for t 0) can show large deviations. We also found that the ratios 

® In fact, we measured the ratios, 


{dMt)Yp{Q)) 


for Q, /3 = 1. • ■ • , 16, 


(4.8) 


In the present notation of the gamma matrices, the ratios (4.8) with a j3 were zero within the statistical 
errors. 
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0.02 



Figure 1. Ratios related to the SUSY WTI at T = 1 for L = 16 and fP = 0.01 (dashed blue line) 
in the SU{2) theory. The horizontal axis is time on which the supercurrents are defined, while the 
source operator is put on t = 0. 


with different spinor indices were well degenerate. In the continuum theory, the sixteen 
supercurrents are interchanged with each other under the transformations. The 

50(9) internal symmetry is broken to SU{2) in the present lattice model. The degeneracies 
indicate that the broken 50(9) symmetry is mostly restored at this lattice spacing. 

We estimated the plateau values using the constant x^-ht, after which we extrapolated 
the fitted values to the continuum limit. Figure 2 shows the extrapolated values for three 
different masses, = 0.01,0.02,0.05. For all of the three masses, the results were about /x^ 
within the statistical errors. Thus, we found that the SUSY WTI holds in the continuum 
limit and the lattice theory surely reproduces the correct continuum theory. 

Finally, we examine the massless limit. For the extrapolation to the massless limit, we 
used the following linear function. 


/(/x^) = c/x^ + d, (4.10) 

where c, d are the fitting parameters. In table 2, c and d are presented for each a. For all of 
the spinor indices, c = 1 and d = 0 within the errors. These results provide the possibility 
that the SUSY theory can be defined by using the extrapolation to the massless limit even 
for the parameter region where the instability is observed. 
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Figure 2. Mass dependence of the plateau values at T = 1 in SU{2). The ar-axis is and the 
y-axis represents the plateau values, extrapolated to the continuum limit, over the mass squared. 


a 

c 

d 

xVdof 

1 

0.997(17)(12) 

-0.00019(32)(17) 

1.7e-04 

2 

0.998(17)(13) 

-0.00022(32)(18) 

4.6e-03 

3 

1.001(17)(15) 

-0.00031(30)(18) 

3.9e-02 

4 

1.000(17)(12) 

-0.00027(32)(21) 

3.0e-02 

5 

0.997(17)(13) 

-0.00021(32)(18) 

3.8e-04 

6 

0.997(17)(12) 

-0.00021(33)(19) 

9.7e-03 

7 

1.001(17)(14) 

-0.00032(30)(19) 

3.9e-02 

8 

1.000(17)(12) 

-0.00028(32)(20) 

3.7e-02 

9 

0.997(17)(12) 

-0.00020(32)(17) 

l.le-05 

10 

0.997(17)(14) 

-0.00021(32)(18) 

6.2e-03 

11 

1.001(17)(14) 

-0.00032(30)(18) 

3.4e-02 

12 

1.000(17)(12) 

-0.00027(32) (20) 

2.9e-02 

13 

0.998(17)(13) 

-0.00022(32)(17) 

9.1e-05 

14 

0.997(17)(13) 

-0.00020(33)(19) 

l.le-02 

15 

1.001(17)(14) 

-0.00033(30)(18) 

4.4e-02 

16 

1.000(17)(12) 

-0.00027(32) (20) 

2.7e-02 


Table 2. Fit results for the plateau values using the linear function, /(^^) = cfi^ + d. 











5 Internal energy of the black hole 


The internal energy of the black hole is a fundamental quantity to study the duality con¬ 
jecture for the present theory. In this section, we see the numerical result of the internal 
energy, where the effect of the phase of the fermion pfafhan is quenched. By comparing 
them with the theoretical prediction of the gravity side, we examine the validity of the 
conjecture. The numerical results of the internal energy are also found in [8], [9], [11], [12], 
[13], [14]. 

The gravity side predicts that the internal energy of the non-extremal black hole, which 
is the gravity dual of the target gauge theory, behaves 

in the large N limit. The leading order coefficient ci can be calculated by hand [5]: 

Cl = ^ (413152 ' = 7.407.... (5.2) 

In contrast, the next-to-leading order coefficient C 2 is unknown. 

In the gauge theory side, one can estimate the internal energy of the gauge theory from 
the expectation value of the boson action [12]. Indeed, for the definition of the internal 
energy (E) = --^logZ with /3 = l/T, we can perform the /3 derivative of logZ analytically 
and obtain 

{E) = -|(S), (5.3) 

where S is given by (2.13). The Q-exactness of the total action implies that the internal 
energy (5.3) vanishes in the zero temperature limit where the effect of the SUSY-breaking 
boundary conditions fade away. 

On the lattice, we define the internal energy using (5.3) with the lattice action (2.17). 
Since the integrations of the fermions and the seven auxiliary fields can be analytically 
performed, it turns out to be 

{E) = |(5s) - ^iVt(iV2 - 1)| , (5.4) 

where k is the number of unintegrated auxiliary fields, and if /c 7 ^ 0 the boson action Sb 
includes the action of the k auxiliary fields. As explained in section 3, the present lattice 
action has the auxiliary field a which remains unintegrated so that k = 1. Using (5.4) with 
the lattice boson action Sb extracted from (2.17) and k = 1, one can measure the internal 
energy from the lattice simulations. 

Figure 3 shows the lattice results of the internal energy, where the red and green points 
denote the results for Y = 14 and ones for N = 32, respectively. The lattice size is L = 16 
for T < 1 and L = 8 for 1 < T < 5 with the corresponding lattice spacing aT^g = 1/L. The 
x-axis denotes dimensionless temperature Tgff = The dashed blue curve represents 
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Figure 3. Internal energy of the black hole against temperature. The simulation results (red for 
N = 14, green for N = 32) coincide with the result of the high temperature expansion (dashed 
orange curve) at high temperature and approach the theoretical prediction (dashed blue curve) as 
the temperature decreases. 


the theoretical prediction of the gravity side at the leading order of the expansion, (5.1) with 
C 2 = 0. The slightly curved orange line is the result of the high temperature expansion 
[26], [27]. As can be seen in the figure, the lattice data coincide with the result of the 
high temperature expansion at high temperature, while as the temperature decreases, they 
smoothly approach the theoretical prediction of the gravity side. 

In figure 4, we focus on low temperature of figure 3. The data surely approach the pre¬ 
diction (dashed blue curve) and are likely to coincide with it as the temperature decreases 
further. But, unfortunately, the temperatures we used in the simulations were not low 
enough to explain the leading behavior of the gravity side. To obtain quantitative results 
for the leading-order term, simulations at further low temperatures are required. 

Instead, one can study the contribution of the next-to-leading order term by fitting 
the lattice results using the following formula, 

f{x) = 7.41x^'® -b Cx^, (5.5) 

where C and p are the fitting parameters. From (5.1), if the duality conjecture is really 
true, the obtained p should be 4.6 within the statistical errors. We performed the fit using 
the 5 points within 0.375 < T < 0.475 and obtained 

C = 9.0 ±2.6, p = 4.74 ±0.35. (5.6) 
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Figure 4. Low temperature region of the internal energy of the black hole. The dashed blue curve 
is again the theoretical prediction of the gravity side up to the leading order. The dashed curve 
denotes the fit result which is obtained by fitting 5 points within T = 0.375 — 0.475. 


The obtained p is consistent with the theoretical prediction of the gravity side within 
about seven percent statistical error. This is the first lattice result of the NLO term, which 
quantitatively shows the validity of the duality conjecture in this system. 

In [9], the NLO term was estimated from the numerical simulation based on the momen¬ 
tum sharp cut-off method using the same fit formula (5.5). The fitted values C = 5.55(7), 
p = 4.58(3) are consistent with our results (5.6) within two sigma. However, those values 
were obtained from their data in a little high temperature region 0.5 < T < 0.7. We 
also tried to fit the lattice results in the same temperature region, but could not obtain 
a reasonable result within x^/dof < 1. This observation raised the possibility that the 
temperature region 0.5 <T< 0.7 was a little high to estimate the next-leading order term. 
The reason of the discrepancy between the two results have not been clear in detail so far. 
Further simulations are now in progress and will give us the final answer. 

6 Summary 

Lattice gauge theory is a promising framework to reveal the gauge/gravity duality and the 
quantum effects of gravity from the gauge theory side. In this paper, we have investigated 
the duality from the lattice simulations of supersymmetric SU{N) Yang-Mills theory in 
l-|-0 dimension with sixteen supercharges. The numerical results of the SUSY WTI have 
shown that the Sugino lattice action that we used reproduces the correct continuum theory 
in the continuum limit. We have also estimated the internal energy of the black hole and 
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found that it tended to approach the prediction of the gravity side at lower temperature, 
and the obtained power of the NLO term indicated the validity of the duality. 

At low temperature with small N, the instability related to the flat directions was 
observed. Whenever we encountered the instability, we changed N to larger ones: = 14 

for Teff = 0.5 ~ 5.0 and N = 32 for = 0.375 ~ 0.5, to avoid the problem. Also, 
all of the results for the internal energy of the black hole are ones obtained within the 
phase-quenched prescription at the fixed lattice sizes; L = 16 for T^g < 1 and L = 8 for 
Teff > 1. Since the lattice spacing aTgff = 1/L was also fixed for each temperature, we have 
to take the continuum limit of the results for the NLO term. For this purpose, the lattice 
simulations at different lattice spacings are ongoing. 
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